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ABSTRACT 

In OTder to reduce the error growth during a nummcal integration, a method of stabilization 
of the differential equations of the Keplerian motion is offered. It is characterized by the 
use of the eccentric anomaly as an independent variable in such a way that the time trans- 
formation is given by a generalized Lagrange formalism. The control terms in the equations 
of motion obtained by this modified Lagrangian give immediately a completely Lyapunov- 
stable set of differential equations. In contrast to other publications, here the equation of 
time integration is modified by a control term which leads to an integral which defined the 
time element for the perturbed Keplerian motion. 

INTRODUCTION 

It is well known that the classical differential equations of the Keplerian motion are 
unstable in the sense of Lyapunov. In general, Lyapunov-unstable differential equations 
develop more unavoidable numerical errors during a numerical integration than Lyapunov- 
stable equations do. We consider here the stabilization of the differential equations of 
Keplerian motion with the aim of improving the accuracy and efficiency during the 
numerical integration. We propose a stabilization method which is purely conservative in 
contrast to other methods (Baumgarte, 1972a). It is characteristic for all conservative 
methods, that they make the revolution time independent of the initial conditions 
(Baumgarte, 1974). 

GENERALIZED LAGRANGIAN 

In this method the stabilization goes hand in hand with the introduction of a new independent 
variable s instead of the time t. This procedure is called time transformation and s is 
called fictitious time. Furthermore, we will require that our stabilized equations of motion 
be developed from the Lagrangian formalism. But, here we have to use an appropriately 
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modified generalized Lagrangian formalism. In order to introduce, instead of t, the 
fictitious time s as a new independent variable, we have to use, instead of the original 
Lagrangian L 


L = L (q., Qj, t). 


I = 1 , 2 ,. . . ,n. 


( 1 ) 


where q. are the coordinates and dot means differentiation with respect to t, the following 
generalized L^rangian L*: 


L* = ht' + p 




( 2 ) 


In equation (2) the prime means differentiation with respect to the independent variables. 
The time t is now a dependent variable (time coordinate). Together with t appears its 
conjugated momentum h, which represents physically the negative total energy. 

Furthermore, we consider only the case, where 

/i = /i(qj,h)>0 (3) 


is a freely chosen positive function only dependent on q. and h. Later we will see that /r 
may be interpreted as the local scale of the time transformation. More general dependences 
of the scale p are also of interest but shall not be considered further here. Our special 
choice of p retains the equivalence between Lagrangian and Hamiltonian formalism. Th' 
choice has two consequences: 

1. The transformed kinetic energy in an inertial system is also a quadratic form in 
the velocity components q' with respect to the fictitious time s. 

2. \ conservative system remains conservative after the time transformation. 


With these restrictions we obtain the differential equations of motion in the following fOTm: 


d 

ds 



dL* 


d 



dL* 


= 0 


= 0 


£ 

ds 





(4a) 

(4b) 

(4c) 


Relative to the system (4a, b, c) the following remarks are in order: 

• Equation (4a) states the equations of motion, which are also given by the 
original Lagrangian L (q., q^, t) after changing from t to s, but which are 
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distinguished only by control terms (Baumgarte, 1972b). This fact will be observed 
li'^er in the special case of the Kepler problem. 

© ' quation (4b) gives: 

h' = p(q.,h) ^ LCq-.ql.h,!), (5) 

dt 

which asserts the energy relation. In the conservative 
ase, where the original Lagrangian L does not depend explicitly 
•n the time t, we get h = constant. 

• From equation (4c) there follows: 

. ( dp d ) 


It will be obvious that the expression in the bracket is a control term, which 
represents the energy relation. This control term vanishes in the exact solution 
of the equations of motion, in such a way that t' = fi. Therefore, it follows as 
previously promised that is the local scale of the time transformation. 


EXAMPLE: KEPLERIAN MOTION 

We consider now the keplerian problem and use Cartesian coordinates. We will choose the 
scale in such a way that the fictitious time s will be for a pure Keplerian motion the 
eccentric anomaly. Only by doing this can we obtain both Lyapunov-stable differential 
equations and equivalence between Lagrange and Hamilton. We, therefore, choose: 


p=— ,r=lxl, (7) 

y/lh 

where r is the ciistanro 

Before we establt.. the generalized Lagrangian L*, we first write the original Lagrangian L. 
With the masa .n = I , K as the gravitational parameter, L has the form: 

I 

L = -lxP+— . (8) 


With the help of equation (2) we now obtain as the generalized Lagrangian L*: 

h 


, y/lh , ,2 

L* = ht + I X I + 

2r " 




(9) 
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With (9) the stabilized equations of motion are: 


(x.5) 



"] 


(lOa) 


h' =0 


(10b) 


5 




Ik^ 

It*' 



(lOc) 


(10a) requires that h* = 0. It is essential that in contrast to an earlier publication (Baumgarte, 
1972b) the vectOT-eq nation for x" does not depend on h. This means the revolution time 
is the fixed number 2rr and this implies Lyapunov-stability for the x" - equation. 

From ( 1 Ob) there follows h = constant. The constant h is computed from the initial con- 
ditions and is placed in the computer as a fixed value. This is the presupposition that now 
in contrast to earlier publications (Baumgarte, 1972b) the time integration (10c) is 
Lyapunov-stable. The expression: 

j_ ( .!£!!• 

2v^l ^ i 

a part of the right hand side of (10c), is proportional to the difference between potential 
and kinetic energy. Therefore, the mean value with reject to the fictitious time s of this 
difference vanishes because we have an oscillator problem in principle. Consequently, the 
expression K} /(2h)^^ in (10c) represents the exact mean value of the right hand side of 
(lOc). Tills fact implies Lyapunov-stability also for the time integration. 


(X)NTROL TERMS AND LYAPUNOV-STABILITY 


In order to show the effect of the stabilization by the control terms, in the classical 
Keplerian equations 


K* 


X = - 




i. 



( 11 ) 


we substitute in place of the independent variable t the fictitious time s by using 
dt/ds = t' = r/y/2h! Because h = constant we obtain: 


(x.x') 


x' + 


2h r 


0 


(12a) 
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h' = 0 


(12b) 


=0. (12c) 

V2h 

We now transpose system (10) in such a way that we can see clearly that system (12) differs 
from system ( 1 0) only by control terms which are the right hand sides of equations (13a) 
and (13c). 


(x-x*) 

X 

Ix'l^ 


1 


X + — 

J.2 ~ 2h 

r 

2r^ 

2hr 

+ — X 
2 ‘ 

(13a) 


h' = 0 




(13b) 

t'- — = 

Ix'l^ 


1 

r 



2r* 

" 2hr 

+ — “ 
2 


(13c) 


In the control terms in ( 1 3a) and (13c) the same bracket appears as a factor, which is 
analytically zero with respect to the energy relation. 

Equation (13c) or (10c) can be integrated. We find as an integral of motion: 



The control terms in (13a, c) produce the Lyapunov-stability under the supposition that 
the constant h is computed once and for all from the initial conditions. The proof for the 
stability of the complete system ( 1 3) or ( 1 0), respectively, (with respect to h = constant, 
whereby this constant is to be computed, finally, by the initial conditions) can easily be 
carried out by making the transformation into action and angle variables, because the 
equivalence between Lagrange and Hamilton exists. By doing this, the generalized 
Hamiltonian, obtained from the generalized Lagrangian by making a Legendre transforma- 
tion, will be linear in the action variables, which implies Lyapunov-stability (Baumgarte, 
1972b). Another proof follows from the equivalence of the system (10a, b) together with 
( 1 4), to the corresponding equations of the KS-transformation (Stiefel and Scheifele, 1971), 

We will call attention to the fact that the dependence of the time transformation 

t' = xfy/^ on the momentum h makes possible the elimination of the instability under the 

presupposition that h' = 0 is integrated exactly. 

In the case of perturbed Keplerian motion, the stabilized system ( 1 3) or ( 1 0) is modified 
by additional perturbation terms. In equation (1 4), C is no longer constant but becomes 
the slowly varying time element (Baumgarte, 1972b; Stiefel and Scheifele, 1971). 
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Numerical experiments have always shown a reduction in error in the numerical integration. 
It appears that the positive effects of the stabilization of the pure Kcplerian motion carry 
over to the perturbed problem. 

We will finally remark that the KS-transformation can be performed directly in the general- 
ized Lagrangiap L* by inserting 




(15) 


thereby giving immediately the analogous complete, stabilized, set of differential equations 
which leads directly to the KS-elements(Stiefel and Scheifele, 1971). 
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